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AN EVALUATION OF CABLE LACING MATERIAL TO SELECT A 
SINGLE TYPE FOR USE ON WIRE BUNDLES 


John P. Tuggle and Leonard G. Rado 
Westinghouse Electric Corporation 
Baltimore, Maryland 


Summary -- A method to assist in the selection of a material for cable lacing 
bundles is discussed. Heretofore little work has been published evaluating cable 
lacings, and opinions vary widely as to the proper material to use. This paper 
is intended to help resolve some of these opinions by presenting the results of 
an experiment performed at our plant. Five rolls of twine and tape were evalua-~ 
ted to determine which material (base material and coating) would best serve our 
requirements. 


TEST PROCEDURE AND ANALYSIS 


In this experiment, five rolls of cord described in Table I were tested for 
breaking strength and elongation before and after sequential fungi-humidity 
tests. In addition, two rolls were tested to determine the effect of a knot upon 
the breaking point and elongation of the cords after undergoing the fungi and 
humidity tests. A diagram of the knot is shown in Fig. 1. Each roll of cord was 
cut into 2) inch lengths, and ten strands from each roll were then selected at 
random for each portion of the tests. These strands were then randomized with 
respect to each other in order to reduce experimental bias. 


TABLE I 
Allowed 
Cross- Breaking Allowed 
Sectional Strength Elongation 
MIL Area (MIL-T-~713) Per Cent 
Roll Material Coating Type Designation (in)2 (pounds) (MIL-T-713) 
ei ed ee ees he eee Pee a A etn bs 
A Nylon Wax Tape MIL-T-713, 0.00096 ite) = 
tape, type (Avg) 
P, waxed, 
Cliae 
B Nylon Rubber Tape om 0.00099 ny = 
C Nylon Rubber Twine ~~ 0.0009 ae =o 
D Nylon Wax Twine MIL-T-713, 0.00080 8 20 
id 4 twine, (Avg) 
Z type P, 
waxed, 
CL 2 
E Dacron Wax. ‘Tape None 0.00096 -- a 


Pie ee 


Upon completion of the fungus test (in which no visible aftereffects were 
noted) the samples were inspected and then placed in the humidity chamber for a 
ten day humidity test per MIL-E-5272. Once again, no visible aftereffects were 
noted. In the compilation of results, a normal distribution was assumed, since 
the data were insufficient to determine the true distribution and since most 
naturally occurring distributions are known to approach the normal distribution 
as the quantity of data (sample size) is increased. 


Figures 2 and 3 show the distribution found for each roll, both before and 
after the fungi-humidity tests. Figure 2 shows the breaking point, while Fig. 3 
shows the elongation. If this test were repeated, the distribution shown on 
these charts would probably vary slightly due to the inherent unrepeatability of 
most test data. However, it can be predicted that the data will lie within spec- 
ified limits 90 per cent of the time. These limits, for the mean and standard 
deviation, are shown in Table II. 


Table III shows the mean tensile strength of each roll before and after the 
fungi-humidity tests. The latter tensile strength is probably inaccurate, since 
the cross-sectional areas of the cords were not remeasured. The tensile strength 
in each case is based on the cross-sectional area of the cord as received. Val- 
ues for roll C are shown to be very high in comparison to the other rolls. Roll 
C is a rubber covered twine and very smoothly finished. Roll D, which is wax 
coated twine, is rather unsymmetrical after the coating process. A check of the 
diameter of the twine before coating (vendor was contacted) revealed that the 
basic twine stranding is the same in both rolls (cross-sectional area prior to 
coating equals 0.000122 square inch). It may thus be assumed that the coatings 
add materially to the thickness and contribute to the very good showing of roll C 
in Tables III and IV. This information must be relegated to secondary importance > 
in light of other results discussed herein. 


Table IV shows the mean modulus at the breaking point for each cord before 
and after test. This modulus will naturally not be the same as Young's modulus 
since the elastic limit of each cord has been exceeded. Also shown in Table IV ~ 
is the ultimate elasticity of each roll. A comparison of data taken before and 
after the fungi-humidity tests reveals a significant difference (99 per cent con- 
fidence or better) between the two with respect to both breaking point and elon- 
gation. The fungus and humidity tests affect the cords detrimentally in both 
cases, except roll A, which is apparently improved by fungi-humidity testing. 
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TABLE IT 


Minimum and Maximum Parameters (90% Confidence Limit) 


Roll Test Min, Mean Max. Mean Min. Std. Dev. Max. Std. Dev. 


Breaking Point (Pounds) 


A Before 59.8 61.2 0.93 2.10 
After 63.0 63.8 0.51 1.15 
B Before 60.0 60.8 0.51 1615 
After 5.0 56.0 1.2; 2.80 
G Before 53.4 54.6 Os 77 1.73 
After 9.7 51.5 1.15 2.60 
D Before 56.2 57.2 0.60 1.35 
After 50.4 55.4 3.14 7-09 
E Before 54.0 54.6 0.35 0.79 
After 51.4 54.8 2.19 4.93 
Elongation (Inches) 
A Before 0.6873 0.7377 0.0316 0.071: 
After 0.6097 0.7279 0.07h1 0.1676 
B Before 0.5376 0.6000 0.0387 0.0883 
After 0.679L 0.7706 0.057) 0.1296 
C Before 0.3967 0.4533 0.0360 0.0806 
After 0.5754 0.6622 0.0538 0.1233 
D Before 0.6777 Oo7i73 0.0435 0.098) 
After 0.7147 0.8229 0.0678 0.1536 
E Before 0.3744 0.4256 0.0319 0.0719 


After 0.4063 0.5155 0.0514 0.1463 


| the experiment indicate that the knot used on samples from rolls 
B ed ee Bn by approximately 2 per cent and reduces elongation by 
approximately 23 per cent. These results apply only to rolls B and D under the 
stated conditions. It may be assumed that the other rolls would be similarly af- 
fected by knots. A diagram of the test equipment used in the experiment is shown 
in Fig. ). Other observations made during the course of the experiment were: 
(1) Twine is more apt to cut a worker's hands than tape} (2) knots in tape are 
easier to untie than those in twine; and (3) the fewer the knots, the longer the 


lacing will last. : 
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TABLE IT1 


Tensile Strength (1b/in) 


Roll Before Test_ After Test 
A 63.02 x (10)2 66.04 x (10)3 
B 61.01 x (10)3 = §5.56 x (10)? 
C 120.20 x (10)3 103.27 x (10) 
D 70.88 x (10)3 66.13 x (10)3 
BE 56.56 x (10)3 55.31 x (10) 


TABLE IV 

Modulus (1b/in®) Ultimate Elasticity (%) 

Roll “Before Test After Test Before Test After Test 
A 8.18 x (10)5 9.13 x (10)9 7.70 7.23 
B 9.92 x (10)? 7.08 x (10)5 6.15 7.81 
C 23.98 x (10)? 15.4 x (10)5 he59 6.69 
D 9.20 x (10)5 7.95 x (10)? 7.70 8.31 
E 13.08 x (10)5 12.10 x (10) 4.32 4.98 


Tepe or Twine Test Specimen 


Drive Motor with 
Graham Speed Reducer 


Pointer Binding Post 
Vi ee ot \ s07>/inse (Dia) 


Output Shaft .75 Inches 
(Dia). Tape or Twine 
winds directly on 
output shaft, 


Vise 


Fig.  - Test setup to measure elongation and tensile pull. 
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CONCLUSION 


If MIL-T-713 is used as a basis for comparison, it may be concluded that all 
specimens are satisfactory, as may be visualized by inserting the specification 
limits in Figs. 2 and 3. Roll A very definitely proved to be the strongest roll 
and is the only roll which remains above the MII-T-713 tensile value before and 
after tensile testing. The superiority of one type of coating over another can- 
not be established from this data. Rolls B and C rate about equal, and roll E, 
which is a wax coated dacron roll, degrades severely in relation to roll A, which 
is a wax coated nylon roll. 


All rolls show significant additional elongation after the humidity-fungus 
tests. Prior to the fungus-humidity testing, dacron proved to be the superior 
material, as borne out by Fig. 3 and the ultimate elasticity figure for roll E in 
Table IV. Based on these findings, wax coated nylon tape is the recommended 
lacing material. The designation for this material will be found in the "MIL 
Designation" column of Table I. 


APPENDIX 


The following statistical methods were used in evaluating the data of the 
experiment. 


Analysis of Variance 


A number of different analyses were made in the course of the evaluation. 
The method used in each case was basically the same. The example given below is 
by far the most complex analysis performed. 


Consider the breaking strength of the cord to be represented by X. Consider 
also that the effect due to environment, vendor, coating, type of cord, and re- 
sidual error are represented by subscripts i, j, k, 1, and m, respectively. A 
given breaking strength with respect to the various effects acting upon the cord 
would, therefore, be represented by Xijxim- If only the effect j were under con- 
sideration, the breaking strength values would be represented by X.4..., or 
simply Xj- Other effects upon the breaking strength would be represented in a 
similar manner. 


The variance s* of each of these effects or combinations of effects upon the 
breaking strength of the cord was calculated as follows: 


Step 1 


5 = isdn Xi jan)* 
nN 


where ™) represents the total number of X values. 


c, ab; skim xy)? 
J 7 jkLm 


where “| jkim represents the number of X values for each effect je 


(i 


where ”] iklm represents the number of X values for each effect j. 


(Gklm X44) 
Cs3 ai eel Eat 9 
i J 7 kim 
where 7) klm represents the number of X values for each group designated by effect 
i and j. 
(film Xjx)2 
Cay =D jk @ilm Xje)* 
J 
"| ilm 


where “) 41m represents the number of X values for each group designated by effect 
peand- Ks 


(&1m Xy Ke? 
Ci5k =yijk Sate 


where 7) Im represents the number of X values for each group designated by effects 
i, j, and k. 


2 
Cyc =Djk1 (sim X41) 
™ im 


where im represents the number of X values for each group designated by effects 
2 
: (em X45«1) 
C3 5k1 Lijkl ym 


where ”) m represents the number of replicate measurements taken under the com- 
bined test conditions. 


Ci jam “Lijkim (X4 5m) 2 
Step 2 
2. Oye =sG 
pica eral 


where p represents the number of different environments. 


2 Cs = 
3 ee 


where q represents the number of different vendors. 


en Nap heme: hee Pent 
S55 x_—e———O 
(p-1) (q-1) 
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Be Cjk ~ Cj 
k(5) Ba (r4-2) 
where r : 
Peek ise baergots the number of measurements in each group designated by ef- 
=e aed eel Ty al let C3 
Me (3) Pg (med) (ped) 


a2 Cjkl - Cik 
ra (jk) x, atte (v, -1) 
peer? represents the number of measurements in each group designated by ef- 
2 S Sajklas Cijk = Cj ~— C5k 
il (jk) ey qr, (v3-1) (p-1) 


Ss Salar eral 
*n CBee weskok ap Pl, Ya-P-a 
The variances ae s 3? s° n a oF and gree ., were compared to the 
: i K(5j)0 i) ik 5) il(kj) 
variance s(43K1) in oe to obtain Snedecor's F value.! ‘The variance Slekray 


was chosen as the denominator on the basis of a method used by C. R. Hicks.2 
The various ratios (F values) were then compared to Merrington and Thompson's 
table3 on the F distribution for significance. 
Comparison and Classification of Means 

Means were compared by a method developed by Tukey. 
Mean and Standard Deviation 

The standard deviations were obtained as follows: 


3 = a x)*/y 
-l 


where s ® standard deviation 
x = measured value 
1 = number of measured values. 


The minimum and maximum means were obtained by 


Misee Sean Nee 
ine v7 


where M ® mean a 
t = a value obtained from Merrington and Thompson's table~ using a signifi- 
cance level of 0.05. 


The minimum and maximim standard deviations were obtained by use of the fol- 


lowing equations: 
s max. ™ S ——— 
Vx (1-1): 0.975 


smin. = S tg 
x" (H-1), 0.025 


where x° is a value obtained from Thompson's table of x2 values.© 
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THE FITTING OF A FAILURE-RESPONSE CURVE TO 
EXPERIMENTAL DATA 


N. M. Peterson 
Convair 
General Dynamics Corporation 
Fort Worth, Texas 


In radiation experiments, effects among components are measured during and/or 
after irradiation with known dosages. These components form a sample from the 
population of components. The fitting of the dose-failure curve is an attempt to 
infer from a given experiment the conditions obtained in a class of components. 
The calculated regression line of the dose-failure diagram is the most accurate 
estimate that can be drawn from the data, granted that the basic assumptions are 
correct. 


To measure susceptibility to radiation, the functional threshold, or 
smallest effective dose, is established. This threshold will vary from sample to 
sample. The dose-failure curve is usually sigmoid in shape; however, a trans- 
formation may be necessary to change the curve to the sigmoid shape. A further 
transformation will establish a linear relation between log of dose and log of 
the ratio of nonfailures to failures, as shown below in Eq. (3). 


If a sample of components of a particular type are irradiated at a constant 
rate and the times recorded at which the individual components cease to function 
properly, the transformed dose-failure curve can be drawn. From these data and 
the regression curve, the functional threshold can be estimated. In a "one-shot" 
experiment, the total sample should be divided into groups, or subsamples. Then 
a separate group of ten or more should be tested at each dose level. The per- 
centage of failures at each dose level will establish the regression line. 


Although the log transform is not an essential part of the treatment, it is 
empirically known to be desirable in many cases. Each new case requires rejusti- 
fication. The final transformed linear relation can be considered as a simple 
case: y = ax + b, where x is log dose. If neither x nor y were subject to any 
uncontrolled, or error, variation, the problem would be easy. But the estimate 
depends on the types of error to which each variable is subject. The variable 
x is subject to errors of dosimetry and y suffers from the many errors which are 
elsewhere considered. 


The precision of the results depends on the variation of the effects 
obtained. To estimate the errors, one must obtain independent estimates of the 
dose used in the experiment and estimate the errors in the parameters involved. 
The response of an individual component may be a continuous variate, or it may be 
an all-or-none relationship. The response of a group of components may be the 
average response or the per cent of components giving the characteristic response. 
The variation in response need not be the same at all levels. 


The characteristics of the dose-response relationship are: 


1. Very low doses elicit no response. 
2. Very high doses elicit total response. 
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3. Between these two doses there exists a point of maximum rate of change 
of response. 


h. The middle range of the dose-response curve is the steepest part of the 
curve. 


DISCUSSION 


If a group of N components are subjected to radiation at a constant rate 
and the failure time for each component is recorded, the per cent of the total 
N that fail in a given time is a measure of the radiation effect. The statis- 
tical distribution of failures in terms of the amount of radiation can then be 


estimated. 


Assume there exists a critical amount of radiation associated with each 
individual component, at which that component fails. Record the time at which 
each component fails while being irradiated at a constant rate of irradiation. 
The components that have failed at any time have a tolerance lower than the 
integrated amount of radiation measured at that time. The parameters of the 
distribution of tolerances can be estimated by comparing the per cent failing at 
each failure time. 


If the total number of failures during the experiment is f, then the total 
fraction failing is f/N, or 100 f/N per cent. If the distribution of failures 
is plotted against the total amount of radiation, the ordinates will. increase 
in uniform steps of 1/N, but the abscissas will fall at irregular points. 


Let x be the total amount of radiation and y the fraction of failures. 
Then, there is a relation between x and y which can be expressed, y = F (x). 
Suppose for simplicity that the function is linear; i.e., y = F (bx ta). If the 
nature of the function can be determined, then the remaining problem is to 

estimate b and a, the parameters of the distribution. 


The principle of linear regression can be applied to this:problem. This 
principle describes the change in y due to arbitrary changes in x. If the plot 
of the failure data is sigmoid in shape, then a transformation is necessary to 
convert the curve into a straight line and simplify the estimation of the param- 
eters of the function. Various methods can be used to determine a transforma- 
tion to accomplish the above purpose. It may be necessary to let u = logio x; 
and then determine y = F(bu + a). 


If y is the fraction of failures at any particular time, then (a - y) is 
the fraction that has not failed. Define the frequency function f(bu + 3) as 
the derivative of F(bu + a); i.e., 


f(bu + a) = F'(bu + a) (1) 
then 
i . 
y= Jes f(bu + a) du. (2) 


If y' = -by (1 - y) where bis a constant, then the solution of the above dif- 
ferential equation is: 
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In (1 - y) -Iny=bu+a 
or 


In==2 = tus+a, (3) 
where a is the arbitrary constant. And if 
eee eee. 
Ms 
then z = bu + a is the linear regression function required. 


Solving this equation for y gives 


1 
VW GouTTaa (1) 


This function is called the "logistics function." It is similar in shape to the 
Gaussian distribution, but simpler in some respects to apply. 


Since the effect of amount of radiation is frequently determined on a loga- 
rithmic scale, it may be necessary to make several transformations. To illus- 
trate: if x is the integrated flux in neutrons/cm2, Jet.t = login. x. A further 
transformation would be to normalize the variables by 

usoet (5) 


where t is the mean of the t's and k is the unit of measure desired. 


To estimate the parameters a and b, the method of least squares can be used. 
The quantity z would be an estimate of a. Since (1 - y) and y are equal at the 
mean, the value of Z is zero in this case. Hence, only the parameter b need be 
estimated. By the method of least squares, an estimate of b is 


Eee é (6) 
Quy 
Another method of estimating the parameter b, which can be done much more rapidly, 
but less accurately, will also be illustrated. The steps in this rapid method 
are: 


1. Justify the use of the logistics function. 
2. Plot the radiation-failure response data. 
3. Select one point near each end of the failure data to use in estimating 


the parameter b. 
lh. Solve for the parameter b, and check solution by comparing with addi- 


tional points on the curve. 


Each of these two methods will be illustrated by an example. The longer, 
more deliberate method should be used where greater precision is desired. 


Method 1 


To illustrate the rapid method of estimation, data from OHNL 1720, Decem— 
ber 6, 1955, are used. These data show the effect of radiation on primers firing 
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by stab action. Figure 1 shows x, the height of fall of a 3.95 ounce ball 
measured in inches as the independent variable, and y, the per cent of primers 


firing. 
The x values can be simplified by the transformation: 
XeeXO 
eens o> 


where x is in inches and u is the number of increments on the abscissa scale. 


The experimental data points for unirradiated primers are tabulated in 
Table I. An arbitrary value of x9 has been selected. The value of u which will 
give 50 per cent response is chosen and designated m. By interpolation 


hte? acts 


The fitted curve will be symmetric about the point m. 


TABLE I 


Unirradiated Primers 


x oot aye 
0.50 ah 1% 
0.75 @) 12% 
Meo) +1 20% 
1.25 +2 80% 
1.50 +3 100% 


To fit a logistics curve to the data in Table I, note that from Eq. (lh) 


fh 


Ai 1+ ebuta 


Aegean eee 22:3 
wg 


Lees EV ee 
y 


Let a = -bm. Then, repeating Eq. (3), 


it =z eb(u-m). 


The value m'= 1.5 can be checked at 
u=1.5 
50 (151, 5) aeons eye a 
“ih 


or 
eyzrl 
y = .50 = 50%, 
which satisfies the condition. To estimate b, choose the point u = 1 
| y= 200% =1/5. 
Then, substituting the value of y in Eq. (3) 
1-1/5 Z 
aos = eb(1 - 3/2) 
= e-1/2 b 
~1/2b = 1.39 
Je) ~ -2.78 
check by using 
Ua= eo 6 = 80% = ,/5 


then 


1 - b/d = i/h = e-2-78(2 - 3/2) = e-1.39 = Wile 


The fitted curve 


za uh 
é Jetbe le 10 Ua~23/o) 


is plotted in Fig. 1. It fits the data closely except at the point x = .75. To 
illustrate a method of obtaining a closer fit, two curves are fitted to the data 
shown in Table II. These data show the effect of irradiation on this type of 
primer. 


TABLE II 


Primers After Irradiation 


Fitted 
x Gib y y 
0.50 =a 0.05 0.05 
0.75 @) 0.52 0.50 
1.00 +1 G77 0.77 
1.25 +2 0.97 0.945 
1.50 +3 1.00 0.973 


Here m is approximately equal to zero. Hence 
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Therefore, if 
uz il 


77 


then 


- .77 


pare Tee 2 3= ede e-1.2 


b = -1.2 


i 
= 2 
Riel rere vo foru#®0. 


A very close fit is achieved by defining the response curve in two sections, 
as shown in Fig. l. 


As aan 


ATT = 


PRIMERS FIRING BY STAB ACTION (%) 


0.25 0.50 0.75 1.00 1.25 1.50 
_ HEIGHT OF FALL OF 3.95 OZ. BALL ON FIRING PIN 
(inches) 
Fig. 1 - Sensitivity change of irradiated M26 primers at ambient 
temperatures. 
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Method 2 


Data to illustrate the least squares method of estimating the parameter b 
are from Admiral tests of radiation effects on semiconductors. Twenty-five 
transistors were tested to determine the fast neutron integrated flux at which 
each transistor failed. There were 22 failures at an integrated fast neutron 
flux less oan 2.7 x 10°7 neutrons/cm*. Let x be the integrated neutron flux in 
Ber oner on si Let t = logjg x and t = average value of t. The transformations 
u= t —- t an 


z= Inst 


lead to the functional relations that are required. The values of x, t, u, 
(1 - y)/y, and z are shown in Table III. Applying the method to estimate b 
gives 
es if 
y at + pool: u “ 


This curve is plotted in Fig. 2 and the failure points are also shown. 


By: plotting the theoretical or fitted distribution of failures curve, it 
is possible to predict the tolerance functional threshold of the group of tran- 
sistors. The functional threshold value below which no transistor would be 
expected to fail is read from Fig. @: 


Rpg x 1019 neutrons/cm@ . 


Further simplifications of the application of Method 2 are possible by plot- 
ting a family of curves on transparent log paper and using these curves as tem- 
plates for estimating the equation to fit the experimental data. 
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Fig. 2 - Transistor failures. 
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TABLE IIl 


Transistor Failures from 25 Transistors Tested 


Failure x = Flux + = 10810 x u = t-t 


Number n/om® 


1 O° OC) * 19301 
2 4.0 0.602 
3 5.5 0.740 
hh 6.) 0.806 
5 7.5 0.875 
6 FeO (12) 12.000 
7 10 0.000 
8 120 0.000 
9 ie 0.079 
10 ieee 0.079 
Vat Ais 0.431 
a 2.9 0.462 
13 3.0 0.477 
1h 3.6 0.556 
aS 4.5 0.653 
16 5.0 0.698 
17 8.0 0.903 
18 9.5 0.977 
19 1.6(13) 13.30) 
20 1.9 0.278 
21 3.0 0.477 
22 hy 


0.568 
Denotes 2 x 1011, 


2 = 273.17141 
= 12.1688 


c+ 


ral 


= 2.611 x 1012 n/em@ 


Hence 
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1d S 
-0.81) 
20 plevie 
=O2610 
-0.5)1 
-0.16 


-0.16 
-0.)16 
-O. 337 
-0. 337 
+0.01h 
+0.0h5 


+0.060 
+0.139 
40.236 
+0.282 
+0.186 
+0.560 


+0. 787 
+0.861 
+1.060 
#1.151 


iy ge ini 
oe ee 
2h 3.178 
2372 2.2 
278 1.992 
21/h 1.658 
20/5 1.386 
19/6 1.153 
18/7 0.910 
17/8 0.753 
16/9 0.571 
15/10 0.05 
14/11 0.239 
13/12 +0.077 
12/13 -0.077 
11/1 -0.239 
10/15 -0.1,05 
9/16 =0.571 
8/17 ~0.753 
7/18 -0.91,0 
6/19 aac153 
6/20 -1. 386 
4/21 -1.658 
3/22 -1.992 


APPLICATION OF A METHOD OF INSPECTION TESTING TO ASSURANCE RELIABILITY 


N. M. Peterson 
Convair 
General Dynamics Corporation 
Fort Worth, Texas 


Summary -- This paper describes the use of a sampling method of inspection test- 
ing for the purpose of controlling reliability in a system. The accep table 
quality level (AQL) and mean failure free time can be estimated by the maximum 
likelihood method, assuming that errors are distributed by the Poisson Law. 
Hence, a choice of quality levels or reliability levels can be made. Single or 
multiple sampling plans may be used depending on the cost factors, the sample 
Size, and the available information on reliability. The sample size, acceptance 
number, and cost due to Type I and Type II errors are each considered. One aim 
of inspection testing is to determine when preventable failures occur. 

An inspection lot in this case is a test run of a system for a given time. The 
system is accepted or rejected on the basis of the number of failures found in 
the sample tests. 


The abject of the inspection system is (1) to accept the reliability of the 
system at least 100 (1-A) per cent of the time when the true life is more than Tj 
and (2) to accept the reliability of the system not more than a small proportion, 
B, of the time when the true mean life is less than To, To < T}. 


DEFINITIONS 


Acceptable quality level. 
Failure-free time, mean life, expected time to failure. 
Estimate of the mean life. 


Acceptance number. A value of time to failure which, if 
exceeded by the sample mean failure free time, will assure 
the required reliability coefficient. 


r Number of failures. 


Inspection Measuring, examining, testing, gaging, or otherwise compar- 
ing the reliability of the system with the required 
standard. 


Type I Error . Vendor's risk. The probability of rejecting the reliabil- 
ity of the system, even though the reliability is 
satisfactory. 


. vai Hy ili ting the relia- 
T Ti.Error Convair's risk. The probability of accepting 1 
ei bility of the system, when the actual reliability is not as 
| good as the specification. 3 | 
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INTRODUCTION 


A weapons system is composed of one or more subsystems. Each subsystem is 
subject to a reliability specification by a vendor and acceptance by Convair. 
It is certain that there will be failures in the system, though the failure rate 
is not known. The user is faced with the problem of setting up standards for 
accepting the vendor's subsystem. He wishes to be assured with a degree of con- 
fidence (1-A) that the reliability of the subsystem will be at least as good 
as R. 


Theoretically, this is the general problem of statistical decision, as for- 
mulated by Wald. The true reliability of the subsystem cannot be determined. 
But there is a random variable T (mean time to failure) whose distribution. 
depends on the true situation. An observed value of T provides information about 
the true condition of the subsystem, and it can be used to increase our chances 
of making a good decision. The general problem is that of choosing a decision 
function d(T), relating the sample space to the decision space. 


Adequate control can be attained by the application of a program of inspec- 
tion and testing. Now, the variety of plans, procedures, and tables for a 
program is unlimited. The conditions affecting a program need to be carefully 
stated, and a plan must be determined which minimizes the maximum risk. This is 
the Minimax principle. This paper presents a plan for testing a hypothesis 
regarding the reliability of a system. A sequential testing program might also 
be devised to determine a system's reliability. 


Before a sampling procedure is set up, several questions should be 
answered: 


1. How long should each test be? 

2. How many subsystems are required to be tested? 

3. What variables will be modified? 

),, What type of sampling should be employed, single or multiple? 
5 


. What is the risk function (the expected loss to the vendor and the user) 
if the Type I error is a function of weight, cost, performance, etc., and 
the Type II error is the probability of catastrophic failure? 


The procedure should determine the degree of the reliability of the system or 
subsystem. The specified reliability of the subsystem should be made binding, 
and provisions should be set forth to be followed to measure such reliability. 


STATEMENT OF THE PROBLEM 


The purpose of this paper is to present a rational scheme for determining 
the reliability of a subsystem as a basis for the decision whether to accept that 
reliability or not. We define inspection as the process of measuring, examin- 
ing, testing, gaging, or otherwise comparing the reliability of the subsystem | 
with the standard required. A failure is any deviation of the functioning of a 
part or subsystem from the requirements of the specifications, drawings, purchase 
descriptions, contract, or order. A critical failure is one that could prevent 
the performance of the tactical function of the system. 
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Assume that the subsystem, before a test run, is checked and serviced, that 
the probability of failure is small, and that failures are distributed by the 
Poisson Law. And consider a random variable, T, the mean time to failure. 

A test run is a trial with the possible outcomes of "failure" or "no failure." 
We can think of this as a Poisson trial since the distribution of failures is 
P(x) = e-U uxX/x! where u = t/T and T = the mean time to failure. The probabil- 
ity of no failures in time t is P(0O) = e-t/T since Of = 1. Hence. the 
probability of one or more failures is Q(0) =1-P(0)=1- e-t/f, and dQ(0) = 
1/T e-t/T dt is the probability of failure in time dt, or the rate of occurence 
of failures. The distribution of a sequence of Poisson trials, given by the 
failure distribution P(x), is thus approximated by the continuous exponential 
function dQ(t)/dt or f(t) = 1/T e-t/T. 


To estimate the mean time to failure T, a sampling scheme is devised. The 
maximum likelihood and minimum chi-square methods of estimation both lead to 
E(x) = X as the estimate of the Poisson parameter u, with minimum variance. 
We might also consider u as a variable of Pearson type III. However, we know 
that x has a Poisson distribution if the parent variable x is distributed by the 
Poisson Law. In the continuous exponential distribution, we wish to estimate the 
parameter T. 


Since we are approximating the Poisson distribution by the given exponential 
function f(t), we investigate the expectation of t: 


co ; 
E(t) sv t f(t)dt = } t/T e-t/T gt s T, the mean life. 
° Ox 


The higher moments can be obtained by using the Laplace transform, since (PK eee 
T (k +1) TK. To verify this equation, let t/T = x, then dt = T dx. Then 
f(t) dt ——» e~* dx and 


CO OD 
{ tk/r e~t/T at — | TK xk e-X dx = TK [“(k + 1) 
(e) [o) 


since 


T(k-1)® 1 xk @-X gx, 


(e) 


The problem then is: given a function f(t,T) = 1/T e-t/T, +> 0, T> 0, deter- 
mine a rule of action to (1) accept the reliability of the system with 
probability 1 - Aif T 2T, and (2) reject the reliability of the system with 
probability1-BifT< 1 where T2 <T . Since we know the moments and the 


distribution of T, we can design this decision function. 


ANALYSIS 


In formating a plan for reliability assurance, the sampling procedure mus t 
be decided, and the risks of either taking unnecessary action or of failing n 
take necessary action for a given shift in reliability mst be determined. he 
acceptance number for a fixed sample size must also be established. 


eu 


A single sampling plan has two parameters, sample size and acceptance num 
ber. We first decide how many systems will be tested and how long each test will 
be; this is the sample size. The failure rate can be determined as a grand mean 
from the total running time and the total number of failures. A reliability 
chart could be designed to show the expected failure rate for "reliable" opera- — 
tion. The primary consideration is the measurement of time between failures. 

The mean time between failures can be easily determined from the times between 
individual failures. We merely divide the total operation hours by the total 
number of failures. 


The sampling scheme can assume a sequential nature, as data are accum- 
lated. However, where in classical theory sample size is determined in advance, 
in sequential analysis the sample values determine the size of the sample needed 
to reach a decision. If a trial run continues to failure, the run length L is 
arbitrary. We can determine the reliability R = e-I/T, the probability of no 
error ty? run of length L. The probability of one or more errors is then 
1 - em ° 


When making a statistical test of a hypothesis, there are two types of 
error. In a Type I error, the test run may have a number of failures greater 
than the limit adopted, although the reliability is satisfactory. The higher the 
limit adopted, the smaller is this probability. Type II error arises when the 
reliability of the system is not adequate, but the number of failures is less 
than the limit set up and no action is taken. This error depends on the shift in 
the value of the parameter. 


Normal inspection is continued as long as the system is functioning in the 
neighborhood of AQL (the average number of failures we expect). If the reliabil- 
ity is consistently better than AQL, the testing program may be shortened. 
However, if failures are numerous, an extended testing period may be necessary to 
give protection against accepting an unreliable product. 


We assume the system is from an infinite population of homogeneous products 
whose failure rate is equal to the process average failure rate which is constant 
during the inspection period. We are willing to accept the system when the true 
failure rate is actually not different from the average number of failures. When 
the failure rate is different, a new test program may be necessary. The average 
failure rate is used as the AQL. We are interested in the mean failure free 
aay the reciprocal of the average failure rate, which is estimated from the 
sample. 


SAMPLE SIZE, ERRORS, AND RISK 


One problem that arises is that of determining the most economical number of 
systems to use in the sample if we are given the probabilities of accepting the 
system, if reliable, and of rejecting it, if unreliable. We can estimate the 
cost of additional testing due to rejecting the reliability of the system. This, 
multiplied by the probability of such an error, gives the expected cost of a Type 
I error. The cost of accepting an unreliable system multiplied by the probabil- 
ity of such an error for a given reliability level gives the expected cost of the 
Type Il error. The sum of these two expected costs gives the total cost for the 


given test program and reliability limits. We would like to find optimum limits 
such that the greatest risk is minimized. 
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Suppose we wish the reliability coefficient R to be 0/70. If R is as good 
as 0.70, we Should like 95 per cent of the time to reach a decision to accept the 
reliability of the system as satisfactory. The Type I error is 0.05, the proba- 
bility of rejecting the product even though it is reliable. If R is less than 
0.70, the reliability is not adequate and 90 per cent of the time we wish to 
reach the decision to reject the product as unsatisfactory. The Type II error is 


0.10, the probability of accepting a system which is not as reliable as we would 
like it to be. 


We can then set up the risk function: risk = 5 per cent of the cost of con- 
tinuing testing plus 10 per cent of the cost of accepting a system which is not 
reliable, plus the cost of the test program. It may be necessary to average the 
expected cost for all possible shifts in reliability level, weighing the average 
according to the frequency of the various shifts. The sum of the two expecta- 
tions gives the total cost due to errors for given sample size and control 
limits. Now we would like to minimize the cost. 


PROCEDURE FOR DETERMINING RELIABILITY 


To determine the reliability of a given system, we want to test the hypothe- 
sis H,: T = Ty, against the alternative Hj: T=T2< TT]. Our rule of action 
is to accept Hp if T > C and reject Hp if T << C. We want 1-A per cent of the 
time to accept H, if true and < B per cent of the time if H) is true. A and B 
are the errors of the first and second kind. It is possible to find values of 
r and C if we are given the ratio Tj/To and A and B. The information necessary 
for making these decisions is given in Table I. 


In many cases we can neglect the first r-l1 failures and just consider the 
time of the rth failure. Therefore, it is possible to set up a sequential sam- 
pling procedure to decide on the length of time that testing must be continued to 
obtain the information necessary to reach a decision regarding the reliability of 
the system. Since 1/T2 - 1/T] > 0, the region of rejection for T = T] is’ 

Tf <C. So we choose C such that Pr (T <C/T=T)) =A. 


We consider the N component parts of the system as a sample of N items drawn 
at random from some infinite population. We test until some component fails. 
The component is replaced by a new component. This is sampling with replacement. 
Each failure will determine the time of the end of a test run. The tests to 
failure can then be ordered, and the first r times to failure can be designated 


as tl, to ° on e tr. 


Epstein and Sobel! have shown that the maximum likelihood estinate of T in 
the distribution f(t) = 1/T e-t/T is 


~ = ty + bles 


isl r 


This estimate is unbiased and has a chi-square distribution which depends only on 
r. The probability density function of f ds 


£8) hE)” rt ort 


T/T or A B C/T1 


3/2 67 0.05 0.05 0.808 
BG 0.05 O10 0.789 
ae 0:10 0.05 0.827 
Fal 0.10 0.10 0.806 

2 Oe 0.05 0.05 0.683) 


19 0.05 0.10 0.655 
18 0.10 0.05 Orr 32 
152" 0.10 3010 0.6866 


5/2 1h 0.05 0.05 0.6046 
sey 0.05 0.10 0.5608 

sia 0.10 0.05 0.6383 

9 0.10 0.10 0.6036 


3 10 0.05 0.05 0.526 
8 0.05 0.10 0.976 
8 0.10 0.05 0.582 
6 0.10 0.10 0.5253 


h [005 05 Coes 
6 0.05 0.10 0.4355 
5 0.10 0.05 0.4865 
h 0.10 0.10 0.4363 

5 Dore Oro r= 0s Goma Ue on 
ewe O50. tO as 0. tuk 
hy 0.10 0.05 0.44363 
3 0.10 0.10 0.3673 


from which we determine limits for the estimates of T. We wish to choose a con- 


stant C such that if T > C we can state that the reliability of the system will 
meet the specification 1-A per cent of the time. 


If we made the transformation 


ae Ty at as T 
Baan? Seon dy 
then f(T) becomes 
1 ort ty Tole orty 1 1 yt _y/2 
ay) = Gar @) = G) Or? Or = 2T(r)e)— @ 


This is the well-known chi-square distribution with 2r degrees of freedom. 
To determine the value of the variable such that P (T< C/T #1) = A, we wish 


P ly Ee 7a Ae 
Ty 


2h 


From the table of the chi-square distribution, we observe the column for the 


1-A level and for 2r degrees of freedom. Set this value equal to 2rC/T} and 
de termine 


To establish the Type II error, the probability of accepting a system which is 
not reliable when it passes the required test, we observe the probability that 
>C/T = To. We wish this error to be less than or equal to B. If 


2 


xX 
T>ce eos 


or. 
then we want to know for what value of To is the probability that 


true, given that T = To. 


Consider P (Y > 2rC/T2) = B. From the table of the chi-square distribution 
for the B level and 2r degrees of freedom, it is found that P (Y > 2rC/To) =B 
holds if 

2 
2 eens tk 
Pras ( te ae 
B 


2 Or 
Then if P (Y > 2rC/T) < B, it follows that 
2 2 2 


TiwtK . 
= 


X 
ete Beer Neg 1A 
2r er Ty % 
B 


for 2r degrees of freedom. See Fig. 1s 


The right-hand side of this expression is monotonic, increasing as r 
increases. We know the left-hand side of the equation is less than unity. Since 
the ratio of chi-squares approaches unity from below, there will be a smallest 
r such that the inequality is satisfied. This will then be the required test as 
outlined above. Table I shows an example of various ratios of T1/T2 with accom- 
panying values of A and B of 0.05 and 0.10, together with the required values of 
r, the number of tests necessary to assure a required reliability coefficient. 


To summarize: we have described the use of sampling methods of quality con- 
trol for the purpose of assuring a given reliability coefficient. Given 
acceptance region A and B and the ratio of T)/T2, it is possible to determine the 
number of tests, R, and the value of GC such that the probability of accepting the 
reliability of the system, if satisfactory, will be at least 1-A and the proba- 
bility of rejecting the reliability of the system, if not satisfactory, will be 


at least 1-B. 
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X° ‘distribution - 
n degrees of 
freedom 


Y; 
Ui ae 


Fig. 1 - Errors of estimating T. 


For example, suppose we wish to establish a reliability of 0.72 for a bomb- 
ing and navigation system which will operate for a ten-hour mission. Then 


=t/T 

R= 

0.72 = enl0/T 

10/Te=—t 3 
Te), 


The mean life of the system must. be at least 30 hours to insure -the required 
reliability. We must estimate the mean life by a testing procedure. The esti- 
mate of T was given previously. 


If we estimate T from our sampling procedure, what confidence can we have in 
the estimate? In-other words, what T is required to assure with a 0.95 confi- 
dence limit that the mean life will be at least 30 hours? Suppose we select as 
an alternative that the mean life is equal to 20 hours. This would imply a reli- 
ability of approximately 0.60. We would then test the hypothesis is ht ae hag 
30 hours against the alternative Hj: T < To = 20 hours. Then, the ratio 
1 / farms 30/20 =. 3/2, 


Consult Table I for A = 0.05, B = v.10, Tj/To = 3/2 for a 0.95 confidence 
limit and for a Type II error of 0.10. There you find r = 55 and C/Ty = 0.789. 
We obtain C = 0.789 x 30 = 24. This gives the desired test information, and we 
determine the following test: accept H, if T > 2h after 55 test runs; accept 
Hy. is 7 < 2). This indicates that the reliability is 0.60 or less and that the 
probability that the reliability is as good as 0.72 is less than 0.10. 
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A RELIABILITY ANALYSIS OF THE EFFECTS OF NUCLEAR RADIATION 
ON THE ELECTRICAL PROPERTIES OF CAPACTTORS 


Re R. Paddock 
Convair 
General Dynamics Corporation 
Fort Worth, Texas 


The mechanism of radiation damage to capacitors is quite complex, and little in- 
formation is presently known regarding the nature of damage from a microscopic 
standpoint. For this reason, radiation effects on capacitors must be based on 
experimental data. Since the field of radiation effects is quite new, available 
data are usually on room temperature or uncontrolled temperature conditions 
during irradiation. Many of the conclusions drawn must then be considered as 
directional only. 


The first atomic explosion in July, 1945, ushered in the dawn of the Atomic 
Age. Since then atomic-powered submarines and power stations have been in suc- 
cessful operation. The greatest challenge, however, is the design of a nuclear- 
powered aircraft. The aircraft designer is now faced with a new environment, a 
new parameter. Nuclear radiation has adverse effects on both electrical and 
mechanical properties of organic and inorganic materials. The purpose of this 
paper is to establish an attitude toward the radiation damage of capacitors where 
mixed reactor fluxes of both gamma photons and high energy neutrons are impinging 
upon the dielectric materials. 


Figure 1 illustrates a typical set of laboratory data of y in terms of x. 
The independent variable x is the irradiation parameter, and the dependent vari- 
able y is the capacitance value at a specific time in the reactor. In other words 
N capacitors are being irradiated and each capacitor has M data points taken. 
Figure 2 shows this data in graphical form. A smooth curve is drawn through the 
data points and is labeled y = f(x). The standard sample deviation is defined as: 


2 
; Hh | (2ye) i So ] : 


This deviation, which is a constant in this analysis, may be used to esti- 
mate a confidence band within the limits of the test data. Referring to confi- 
dence tables which may be found in most mathematics handbooks, a value 1 is 
chosen which is a function of probability and sample size. The value of 1 sets 
an upper and lower limit on the confidence band. The observer is ¥Y per cent con- 
fident that the variability of the data lies within this envelope. Since most 
capacitors have tolerance limits, a maximum and minimum limit is placed on the 
data. In this case the capacitance exhibits a general decreasing trend so only 
the minimum limit is employed. As the edges of the envelope exceed this limit 
a pessimistic threshold (T,) and an optimistic threshold (To) are established. 
The pessimistic threshold is the maxinum amount of irradiation the capacitors 
will withstand at Y per cent confidence to yield a reliability of Re The op- 
timistic threshold Tp will not be considered since it has no intrinsic value. 


oe 
een 


In Fig. 3, as different values of reliability are chosen, the distance Is 
varies and intercepts the engineering limit along the abscissa. Figure ) shows 
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Fig. 2 


these values of reliability plotted against the capacitor's functional threshold. 
It is interesting to note that the functional threshold curve is similar in shape 
to that of the raw data points. High quality capacitors have been found to have 

a low variability among the individual data points yielding high values of relia- 
bility during irradiation. 


It should be emphasized that the analysis technique does not attempt to pro- 
vide an exact appraisal of the reliability of a component; rather it provides a 
certain degree of confidence that the reliability of a capacitor will not be be- 
low the calculated curve. As such, the reliability shown here is not only a 
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ENGINEERING 
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ENERGY ABSORBED IN DIELECTRIC 
(ergs/gm) 


Fie se3 


function of the performance of the component but also a function of the number of 
components tested; an increase in the number of specimens would raise the relia- 
bility curve. 


A large variety of materials are employed in the construction of capacitors. 
Both organic and inorganic materials are used as the dielectric. The plates of 
capacitors are metallic and are generally radiation-resistant with respect to 
other materials. The integrity of capacitors is in large measure dependent upon 
the vulnerability of the dielectric. 


J % CONFIDENCE 


% RELIABILITY 


FUNCTIONAL THRESHOLD 
(ergs/gm) 


Figs 2 
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Functional threshold determination must recognize both physical and elec- 
trical degradation as caused by irradiation. For example, the electrical prop- 
erties of a capacitor may remain within design specifications but fail because 
of structural degradation in a nuclear environment. Factors which affect di- 
electric materials in nuclear environments are accelerated aging due to changes 
in conductivity. High energy radiation produces ionization in the dielectric, 
flashover, lower breakdown voltages, and lower values of leakage resistance. 


The material is internally heated by energy absorption from the radiation 
field. This energy absorption may be calculated from the following equation: 


Ey * (to Cy + dy Cy)t 


where E, is the energy absorbed in ergs/gm 
g is the radiation flux in gammas or neutrons/cm?-sec 
C is the energy absorption coefficient in ergs/gm 


gammas or neu trons/ cme/ sec 


t is the time in seconds. 


The energy absorption coefficient which was developed at Convair by Brown and 
Fink is the subject for a paper itself, so we will define it as a conversion fac- 
tor of irradiation flux to energy absorbed in a particular material. 


Figure 6 illustrates the nuclear and other irradiations more commonly en- 
countered. The neutrons which are present in the environment have velocities 
which range from slow to intermediate to. fast. You will notice from the numbers 
shown here that the corresponding range covered in energy is tremendous, ranging 
from essentially zero kinetic energy to 10 million electron volts. The electron 
volt is a widely used unit of energy and is equal in magnitude to the gain in 
kinetic energy of an electron falling through a potential of one volt. Gamma ir- 
radiation, which is electromagnetic in nature possesses quantum energies ranging 
from 103 to 107 electron volts. X-radiation and visible light are forms of elec-- 
tromagnetic radiation familiar to all of us. From this summary it can be seen 
that they differ from gammas in that they have much lower energy levels. 


A brief summary of some of the terms we use in discussing this environment 
is shown in Fig. 7. We frequently speak of the flux of nuclear radiation -- 


RADIATION EFFECTS IN DIELECTRICS TYPES OF RADIATION 
OF CAPACITORS 
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RADIATION UNITS 


Flux --------- (particles or energy) /cm2-sec 


Integrated flux -- (particles or energy) /cm2@ 
Dose rate --------~.~---.~-~-~---~-~-~ 


com Ge ED © Se om om ee es OD a OO ee ee ee oe em 


Energy absorbed ------~------ 


Energy absorption rate ----------- Rad/hr or 
ergs/gm-hr 
AEC tolerance dose rate ---~----- 0.3 Rem/wk 


Fig. 7 


either gammas or neutrons -- where we are referring either to the number of nu- 
aoe particles/cm¢-sec impinging upon a receiver or to the amount of energy/ 
cm‘=-sec impinging upon these receivers. If a flux is integrated with respect to 
time, the quantity referred to as the integrated flux is obtained, which simply 
represents the total number of particles/cm2 or energy/cm¢ which is incident upon 
the receiver. From these considerations we say that our receiver, whether human 
or inanimate, is receiving a certain dose of irradiation. 


We measure the rate at which this dose is being given in a unit known as 
rem/hr, and a dose received as a result of prolonged exposure we measure in rem. 
This dose, of course, is simply the dose rate multiplied by the time. Not all 
the energy which falls upon a receiver is absorbed. To measure that energy which 
is absorbed, another unit is used, namely the rad. By definition the rad is 100 
ergs of energy absorbed per gram of material. In order to give a feeling for the 
size of these units, the Atomic Energy Commission has decided that personnel 
working in radiation fields for lO hours/week can receive 7.5 millirem/hour with 
no ill effects. x 


Deleterious effects can be divided roughly into two categories. Some occur 
immediately when the component is placed in a nuclear environment and are sensi- 
tive functions of nuclear irradiation flux density; these are termed rate ef- 
fects. Long term degenerative effects associated with the total irradiation are 
called integral effects. The damage to a given component will depend upon the 
modus operandi during irradiation either static or dynamic, the materials com- 
prising the component, the intensity and length of the radiation period, and 
other environmental factors such as reactor temperature and humidity. 


Figure 8 shows the leakage resistance of several capacitors as a function of 
the energy absorption rate in the dielectric. Vendor and military specifications 
set a minimum leakage resistance for each capacitor. When this leakage resistance 
exceeds this value, the aircraft designer consicers it a failure. A reliability 
analysis was not used here since only one capacitor of each type was irradiated. 


Consider a plan view of a nuclear-powered aircraft with isodose-rate lines 
due to the influence of the reactor for both gamma and neutron fluxes. Assume 
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that a capacitor in question to be placed in a particular electronics package has 
previous radiation curves which are available to the aircraft designer. He may 
predict whether or not this capacitor will exceed its specification limit due to 
rate effects of the flux or to integral effects over a specified flight period. 
From previous discussion we may conclude that the leakage resistance of the cap- 
acitor is dose-rate sensitive and the capacitance is integral-rate or dose sensi- 
tive. It is, therefore, obvious that considerable data must be available encom- 
passing each type component to be used in a nuclear-powered aircraft so that the 
designer may judiciously position the various systems in the aircraft. 


LEAKAGE RESISTANCE (kilomegN) 


ENERGY ABSORPTION RATE IN DIELECTRIC 
(ergs/gm-hr) 


Fig. 8 


The frontiers of space no longer live in the realm of dreams. Eventually 
aircraft guided and propelled by nuclear energy will leave the earth to venture 
across interplanetary distances. From the hearts of atoms the energy of the 
universe has been released on earth. But whether nuclear energy will eventually 
destroy our present world or help to create a new and better one lies not in the 
hearts of atoms but in the hearts of men. Human hate and ignorance or human love 
and knowledge are the masters. The atom will survive either. 


} ACKNOWLEDGMENT 
The data and preparation of this’ paper were accompli 
nd plished while in the emplo 
of Convair, a division of General Dynamics Corporation, in the Nuclear ingineeioat 


ing Department, Design Studies Group. This employment was in. con, ; j, 
: wet unction with 
the cooperative student plan of Southern Methodist University. junction wi 


32 


The aid of W. F. Ekern, J. H. Levine and Hilda H. Hooper in the preparation 


of text and data, and Mrs. J. P. Beck in typing the manuscript is gratefully 


acknowledged. 

BIBLIOGRAPHY 

1. Copp, M. P.: "The Atomic Revolution." General Dynamics Corporation, 1957. 

2. Paddock, R. Re: "The Effects of Nuclear Radiation on the Electrical Proper- 
ties of Capacitors." Convair, Fort Worth, Nov. 17, 1957. 

3. Secrest, E. L.: "The Effects of Nuclear Radiation on Electronic Components." 
Report FZM-701, Convair, Fort Worth, Feb. 19, 1957. 

he. "How radiation effects materials." Nucleonics (Sept. 1956). 

5. Fink, W. L.: "A Method for Comparing Radiation Effects Data From Various 
Facilities." Report FZM-1013, Convair, Fort Worth, 1957. 

6. Knowlton, A. E.: "Standard Handbook for Electrical Engineers." New York, 
McGraw-Hill Book Co., 8th ed., 1949, pp. 522-535. 

7. "Researchers detail progress on many fronts in effort to solve environment 
problems." Military Electronics, vol. 3, no. 1 (July 1957). 

8, Westman, H. P.: "Reference Data for Radio Engineers." New York, American 
Book-Stratford Press, lth ed., 1956, pp. 83-106. 

9, "Solid State Division Semi-Annual Progress Report for Period Ending Febru- 
ary 28, 1955." Report ORNL 1852, Oak Ridge National Laboratories, June 1955 
(Unclassified). 

10, Fink, W. L.: "Methods and Procedures for Radiation Damage Analysis." Report 
MR-N-169, Convair, Fort Worth, June 1957 (Unclassified). 

11. Systems Panel Test #2, Raw data taken at Convair, Fort Worth, Nov. 1956. 

12. Praff, E. Ror "The Effects of Nuclear Radiation on Electronic Components." 
Report #4, Admiral Corp, July 1956 (Unclassified). 

13. Collins, C. G., and Calkins, V. P.: "Radiation Damage to Elastomers, Organic 
Fluids and Plastics." Report APEX-261, General Electric Atomic Products 
Division, Sept. 1956, 13, 66 (Unclassified). 

ly. Pfaff, E. Ro: "The Effects of Nuclear Radiation on Electronic Components." 
Report #5, Admiral Corp., Oct. 1956 (Unclassified). 

15. Pfaff, E. Ret "The Effects of Nuclear Radiation on Electronic Components." 
Report #6, Admiral Corp., Jan. 1957 (Unclassified). 

16, Miglicco, P. S.: "Effects of Reactor Radiation on the Electrical Properties 


of Electronic Components." Report MR-N-116, Convair, Fort Worth, Dec. 1956 
(Classified). 


25 


MULTIENVIRONMENTAL LIFE TESTING OF PARTS AND COMPONENTS 
IN ROCKETS AND GUIDED MISSILES BY STATISTICAL DESIGN 


D. Kevan Lloyd 
The Ramo-Wooldridge Corporation 
Los Angeles, California 


Summary -- This paper suggests a method of studying the effect of environmental 
factors on missile parts and components for the following cases: 


1. When the effects of preoperational environments are being investigated. 


2. When environmental test facilities are sufficiently sophisticated to allow 
simultaneous application of operational environments to be made while the 
item is functioning. 


3. When the first two conditions are satisfied and both types of information 
are being studied. 


The investigation of the behavior pattern of the items over the environmental 
space is carried out by combining censored life testing with statistical experi- 
mental designs and illustrating the results by response surfaces. 


THE GENERAL PROBLEM 


It is well known that because of the complexity of a missile, any part or 
component, etc., associated with it must necessarity have an extremely high reli- 
ability number. To demonstrate or investigate the reliability of a component is 
very difficult. To show that the component will function for the minimum time 
without malfunction requires such a large sample size that it is impractical. 
Even if this could be done, the reliability demonstrated would be true only for 
the conditions and environmental combination under which the test was run. In 
addition, a large number of environmental combinations are possible. Higher 
level system interactions further complicate and vitiate the value of component 
reliability estimates. Finally, because of the number of components involved, 
each must have an extremely high reliability target. 


At the part level, and even at the component level, it is not particularly 
realistic to attempt to demonstrate an absolute reliability number because the 
subsystem interactions cannot be easily generated or simulated when testing iso- 
lated components. It is more satisfactory to subject the components to a variety 
of environmental factors so that the best of several types of components and 
their reactions to the environment can be determined, Component testing may be 
regarded as a means of improving reliability; subsystem testing, a way of deter- 
mining a realistic prediction of reliability and guiding component improvement by 
systematically discovering component problems and weaknesses within the higher 
systems and thereby establishing a natural priority of reliability apportionment. 
However, when it is established that a component or part is critical within a 
system, it then becomes necessary to discover the cause of the trouble by inves- 
tigating the environmental factors and the intensities of the factors which cre- 
ate this situation. If it can be shown that a Significant correlation exists be- 
tween failure rates and certain environments, then appropriate corrective action 
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can be taken. It is the purpose of this paper to suggest how the correlations 
may be established. 


LIFE TESTING AND EXPERIMENTAL DESIGN 


Life Testing 


Life testing, or test to failure, reduces the sample size by extracting a 
maximum amount of information from each sample and instead of measuring an at- 
tribute (success or failure) we are measuring a variate, time to failure. Var- 
iate sampling, being more sensitive and informative, is preferable, Life testing 
demonstrates how long the part or component will withstand the environmental 
stress compared with the length of operation required of it. However, this is 
normally done for one environment at a time, which is essentially the classical 
method of experimentation; i.e., it does not take into effect the influence of 
the environments which, in operational use, act and change their groupings and 
levels simultaneously. Therefore, in order to arrive at a realistic idea of how 
the different factors and their different levels of intensity influence the be- 
havior of the items being tested, the techniques of statistical experimental de- 
Sign and analysis should be applied. 


Because of test facility availability, expense, etc., it is not always de- 
sirable to continue each test until the unit fails; thus the slightly more con- 
plicated method of censored life testing will be discussed in this paper. The 
methods suggested, however, apply equally well to "simple" life testing wherein 
all units are tested to failure. Censored life testing is the name given to 
the method of testing for a given length of time and then terminating the opera- 
tion, having measured the length of operation of each of those units which failed 
before "shutdown" and the number of units which were still functioning at time of 
termination. 


The estimation of the life distribution is carried out as follows: Consider 
the estimation of the life distribution under the environmental combination 


(Az Bj Cisco) multike..)'s 


It is assumed that the life distribution is described by the exponential distri- 
bution and is given by 


f(t) = Wijk... CxP (-Hijk... 4) (1) 


where p33, is the mean failure rate for environmental combination (ijk...). 
Test N eaten instead of allowing all items in the sample to fail, terminate 
the test after time T. The problem is to estimate pjijx,,, - Define the 
following: 


T = termination time 

N = total number of items in the sample 

n = number of items failing before time T 
t, = time of failure of the rth unit (r= n). 


It can be shown! (see Appendix A) that the maximum likelihood estimator of 
p is given by ; 
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n 


n 
ty + (N-n) T 


and x pe ( 2 ) 
8 feel po Be) rari a 


N(1 - e-b?) 


Therefore, under any environmental combination applied to N items for time T, by 

observing n and so ty, and computing Hijk... from Eq. (2), the life frequency 
r=o 

distribution for the environmental (ijk...) is completely defined. 


Experimental Designs 


Now suppose a censored life test be carried out for various values of i,j,k, 
eee, then for each combination we can estimate a life distribution in the manner 
described above. However, if the values of i,j,k,... are properly chosen; i.@., 
the correct combinations are determined by the considerations followed when es- 
tablishing a statistical design, then a table of pijk.., for all relevent values 
of i,j,k,+.. can be formed which will be similar to Table I and which can be used 
in two ways. 


Firstly, a response surface can be generated, thus illustrating the behavior 
characteristics of the tested item over the environmental space; and secondly, 
after a suitable transformation2 (see Appendix B) of the data within Table I, the 
techniques of analysis of variance can be utilized in order to determine signifi- 
cance levels and allow confidence surfaces to be plotted. The testing, there- 
fore ,becomes a statistical design replicated N times with the same termination 
time T for each combination. In order to illustrate the foregoing statements and 
also prevent the notation from becoming too cumbersome, a three-factor environ- 
mental example will be given. 


EXAMPLE: THREE-FACTOR ENVIRONMENTAL STUDY 


A study is to be made of the effect of the preoperational environmental fac- 
tors upon:a component. The environments considered in this study are vibration, 
humidity, and temperature cycling -- factors which would normally be experienced 
by a weapon system transported to a tropical area and stored prior to operational 
use. The following environmental levels were deemed to be representative of the 
expected operational stress: 


Vibration a a] ao 
Humidity bo by 


Temperature Co Cy. 
cycling 


There are 3 x 2% = 12 combinations of environmental factors and levels. If the 
study were carried out as a simple factorial design replicated N times, then, 
having subjected each group of size N to one particular environmentel combination 
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as dictated by the statistical design, all 12 groups of components would be oper= 
ated for time T, and the numbers of components failing and the times at which 
they failed noted. For each group or environmental combination (ijk), (1=0,1,2; 
j = 0,13; k = 0,1), we can compute 


M4 5k 


EN ty + (N - nyj_) T 


1g Le) 


where N and T are common to each group. We can then form the array given in 
Table I. 


TABLE I 


Mean Life Estimates Under Applied 
Environmental Combinations 


Vibration 
ae Sate Ua a 
ature Humidity 
Gate Ged 5 Tat bon ahi 5 Pov tes 
co Pooo #010 F100 F110 #200 F210 
Cy Pool Holl #101 #112 201 #211 


Response Surface Example 


The parameter p44, determines the life distribution of the component under 
environmental combination (ijk). Therefore, from Table I we can generate envi- 
ronmental time-to-failure response surfaces similar to the one shown in Fig. 1. 
This diagram is constructed as follows: take the mean of the p's for each level 
of vibration, which gives the expected value of p44, for ap as 


22 Pojk 
rape a sy 
Po. e hy 
2 Hijk 2 Pagk 
Similarly, for a, and ap the expected values are ass an pee which are 


Hl.. and po, .» respectively. 


Now these values of p for each level of vibration are averaged values over 


h life curve we 
he remainder of the environmental space. Consequently, for eac 
ee a more realistic sample of the behavior characteristics of the components. 
This differs from the classical or traditional method of life testing where only 
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the environmental factor being studied has its levels changed. The construction 
of a frequency environmental response surface can be carried out for each factor, 
graphically indicating whether there is any correlation between failure rate and 
environmental level. 


Probability of Successful Operation 
If top is the minimum required operating time for the component, the probab- 


jlity of it operating longer than tp , the component having experienced a vibra- 
tional level of a, is given by 


eS sitios a Hojk 
P (ap,to) = 2 7 exp - ca t | dt 
to 
FHosk 
= exp |- ‘ to J e 


The surface generated and shown in Fig. 1 is a "template frequency surface" 
because the levels of the factors are discrete. However, if sufficient levels 
are studied and a regression equation between the failure rate and environmental 
factor can be determined, p = g (a), say, where a is vibration and p is the fail- 


FAILURE RATE 


MINIMUM PERMISSIBLE 


OPERATING TIME OF -» 7 
COMPONENT 


Fig. 1 - Vibration frequency response surface. (The shaded area 
represents the expected failure region. ) 
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ure rate, the probability that the component will operate for a time greater than 
to, having received a vibration intensity of level ay, is given by 


P (ayy to) = exp ic g (ay) to |. 


This, then is the probability of successful operation having received a level a, 
of vibration and an average effect of all the other factors studied. 


Establishing the Relative Importance of the Effect of the Environmental Factors 
ce ee ee Eee OL We nV ronmental Pacuurs 


The above method of presentation is effective in illustrating the actual and 
potential trouble region of the components in the environmental space. However, 
much of the value of the method would be lost if it were not possible to derive a 
measure of variability so that significant differences between factors and their 
levels can be determined. It is for this reason that the environmental combina- 
tions were chosen so as to form a statistical designs i.e., an array of the type 
illustrated in Table I. However, the data as it stands in Table I is not in a 
form suitable for analysis of variance techniques, and a transformation is nec- 
essary (see Appendix B) to vitiate the relationship between p and var ee 


After the transformation has been performed the analysis of variance tech- 
niques can be applied. The significance of the effects of the different environ- 
ments and their levels on the lives of the components is measurable, and we now 
have a quantitative basis for directing corrective action. Interactions between 
two or more environments which cannot be detected using the classical method of 
analysis will be indicated by the analysis. The example used for the statistical 
design is a simple factorial; however, the type of design used will depend on: 


1. the number of environmental factors to be investigated 


2. the number of test levels of each factor and the range of intensity of 
each factor 


3. whether it is known or anticipated that there are important interactions 
between some of the factors concerned 


lh. the number of units available for testing and cost of testing 


j : i ivi i d the very 
Much ingenuity can be exercised in arriving at the optimum design, an 
erin tools of the analysis of variance can then be utilized to understand the 
life behavior characteristics of the component over the environmental space. 


A Practical Consideration When Applying the Technique 
ce ee LN ee 


i discussion is concerned with the investigation of the item's 
fens See nee over a range of environmental conditions. However, 
it may not be necessary to fully investigate the behavior characteristics af dur- 
ing the course of the experimenting it becomes apparent that the item is reliable. 
For instance, we are interested in finding the areas of failure and therefore we 
want some of the items to fail; consequently, it is suggested that the most eX= 
treme of environmental combinations (i.e., agbjc, in the example) be mere : 
first. The reason for this is that if no units fail, or very few fail, indicat- 
ing a low failure rate under the most extreme conditions, then by the theory of 
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peripheral testing we may extrapolate back to the nonextreme conditions and 
reasonably assume that the failure rate would be smaller. For instance, if we 
could state with a high degree of confidence that 99 per cent of the items could 
be expected to exceed the minimum operating time under the extreme condition, 
then this might be a sufficient guarantee that the item was satisfactory for all 
conditions without further testing. When the item or testing is expensive, it 
would not be prudent to spend more time or money than is necessary to assure our- 
selves of the item's competence. 


However, if no firm decision can be deduced from the testing under extreme 
conditions, the "demonstration" necessarily becomes an investigation of the 
item's behavior characteristics, because the theory of peripheral testing can no 
longer be applied. Consequently, the sampling should progress from the extreme 
conditions in such a manner that the area of the environmental space under inves- 
tigation is "self-sufficient" for a statistical analysis. Figure 2 and the sub- 
sequent explanation wil] clarify the foregoing statements. 


Consider a two-environmental example for purposes of illustration. The most 
severe combination of environments is ab). If reliability cannot be inferred 
for the remainder of the space from the N items tested under conditions a3b),, ex- 
tend the area of investigation to environmental combinations a3b3, a9b3, agb), and 
perform the statistical operations as described in the main body of the paper. 
This operation is extended as necessary over the environmental space until suffi- 
cient information has been derived to arrive at conclusions concerning the ade- 
quacy of the reliability, the effect of environments, and the possible need for 
design changes. 
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Fig. 2 - Progression of testing. 
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APPENDIX A 


The frequency and probability functions of the life distribution are assumed 
to be exponentials i.e., 


f(t) = pe7Pt, O<tsT, 
P=P (t>T) se Fl pret<o, 


where T is the termination time of the censored test. Taking logarithms and ob- 
taining the first and second differentials with respect to hn, 


OJloreteest ae log f . -1 
log f = log ETE ete See ens Peay apap ember = 2 
Daa 31 All mi ’ acer ne 
@ log P 8° log P 
log P = - nT, Be ented Weed She yy Sd Sav aad eS 8) 
Ou one 


NC, p® e7P [st + (N-n) T), noo, 
The likelihood function L # 
e~PNT > ni =.0, 


where N is the number in sample and n is the number failing before. time T. 


The likelihood estimator of p is given by: 


giving " = 


The likelihood estimator variance p is given by: 


oO 
1 8? log g(@)\ 2(e)de 
var 5 ee oc teelee ns 

= an @=9 


where g(6) is the frequency function of the variate 6.3 In this case, the equiv- 
alent expression is: 
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APPENDIX B 


In order to bring Table I into a form suitable for analysis of variance 
techniques, it is necessary to transform the data so that the means are inde- 
pendent of the variances.® The transformation is given by 


m 


No exact solution was found for the indefinite integral shown, and the ap- 
proximate solution is dependent on the value of pT. When pT is large the .trans- 
formation is logarithmic; when pT is small the square root transformation is ap- 
propriate. However, these two transformations are useful in indicating the type 
of transformation required. The correct transformation will be determined by the 
data generated by the experiment. It is sufficient to find only an approximate 
transformation since the analysis of variance loses little power until the data 
deviate appreciably from normality. 
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